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Abstract. Thermodynamic properties of the spinless Falicov-Kimball model are studied on a triangular 
lattice using numerical diagonalization technique with Monte-Carlo simulation algorithm. Discontinuous 
metal-insulator transition is observed at finite temperature. Unlike the case of square lattice, here we ob- 
serve that the finite temperature effect is not able to smear out the discontinuous metal-insulator transition 
seen in the ground state. Calculation of specific heat (Cv) shows single and double peak structures for dif- 
ferent values of parameters like on-site correlation strength ([/), /—electron energy (Ef) and temperature. 

PACS. 71.45.Lr Charge-density-wave systems - 71.30.-|-h Metal-insulator transitions and other electronic 
transitions - 64.75.Gh Phase separation and segregation in model systems 
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1 Introduction 

The Falicov-Kimball model (FKM) [T] is one of the sim- 
plest and most successful model for the correlated elec- 
tron systems. This model was introduced to describe the 
semiconductor-metal transition in mixed- valence compou- 
nds and rare-earth systems. Recently systems such as tran- 
sition-metal dichalcogenides [1113111] , cobaltates , Gdl2 O 
[7] and its doped variant Gdl2Hx [8] have attracted con- 
siderable attention as they exhibit remarkable cooperative 
phenomena, like metal- insulator transition, charge and 
magnetic order, excitonic instability [2] and possible non- 
Fermi liquid states [Tlin]- These systems are characterized 
by the presence of localized and itinerant electrons con- 
fined to the two-dimensional triangular lattice. 

In the FKM effective interactions are mediated by band 
electrons [TUlfTT] . Therefore, the model has also been stud- 
ied as a model of crystallization, for example, as a model of 
binary alloy. The model describes the systems exhibiting 
different types of phase configurations e.g. regular phas- 
e [H] , phase separation [T^fHlfTB II 1 6 lITT] and stripe phases 
I191I20| . There are predominantly more rigorous results avail- 
able for the FKM on a bipartite lattice than on non- 
bipartite lattices. On a bipartite lattice, one of the most 
important result proved by Kennedy and Lieb [TOl[TT| is 
that at low enough temperature the half-filled FKM pos- 
sesses a long range order, i.e., the localized electrons form 
a checkerboard pattern, the same as in the ground state. 
This result holds for arbitrary bipartite lattices in dimen- 
sions d > 2 and for all values of on-site Coulomb correla- 
tion strength U. 

Recently few interesting results are reported for the 
phase diagram of localized /—electrons due to the inclu- 
sion of correlated hopping (f) [HllllllllMllSSllSnillZ] of 



itinerant electrons in FKM. Phase segregation, even in 
weak correlation limit {U ^ 0.5 or 1), is one of the most 
striking outcome of the generalization of the FKM on bi- 
partite P^[^[^ and on non-bipartite lattices. 

Even though the FKM is one of the simplest models 
for correlated electrons, a clear picture of metal-insulator 
transition in FKM is still debated for many years. Metal- 
insulator transition strongly depends upon the choice of 
the d-electron density of states (DOS). There are only few 
exact results available on metal-insulator transition on a 
non-bipartite lattice [27.29 , 30^ . In a recent work [27| we 
have studied the ground state properties of the Falicov- 
Kimball Model on a triangular lattice with correlated hop- 
ping and observed first order phase transition in the f{d) — 
electron occupation at particular value of /—electron en- 
ergy Ef. This shows that the small change in the lattice 
parameter by means of introducing impurity or by apply- 
ing pressure one can successfully explain the valence and 
metal - insulator transitions. In a separate work |30j we 
studied the effects of Coulomb repulsion between f- and d- 
lectrons (U) and also between f-electrons(C//) themselves 
in a two-fold degenerate f-level on the ground state prop- 
erties. We observed first order (discontinuous) insulator- 
metal transition at particular values of U and Uf. 

A very important question, in this context arises as 
to what is the nature of these phase-transitions at finite 
temperature? Does this model again give the first order 
(discontinuous) insulator to metal transitions of f(d)- elec- 
tron occupation with Ef and temperature? In addition it 
is also interesting to know the dependence of the specific 
heat (Cv) on temperature for different range of parameter 
values. 

In this work, therefore, we study the FKM in its spin- 
less version numerically for all ranges of interactions, and 
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Fig. 1. (n/ - Ef) phase diagram for different values of temperature at U=Q and temperature dependence of the f-electron 
occupation number n/ at (7 = for different values of Ef (shown in inset). 



explore the metal-insulator transition and specific heat on 
a triangular lattice at finite temperature. 

The Hamiltonian of the system may be written as 

+EfY.flf^ + UY,fU^4d,. (1) 

i i 

Here d\,di{fl , fi) are, respectively, the creation and an- 
nihilation operators for itinerant d-electrons (localized /- 
electrons) at the site i. The first term in Eq.(l) is a mea- 
sure of kinetic energy of d-electrons on a triangular lattice: 
only nearest-neighbor hopping is considered. The second 
term represents the dispersionless energy level Efoi the 
/-electrons while the third term is the on-site Coulomb 
repulsion between d- and /-electrons. 



2 Methodology 

In the Hamiltonian (1), hybridization between d- and f- 
electrons is absent. Therefore, local /-electron occupation 
number fij^i ~ fifi, is conserved. The uji = f^fi is a 
good quantum number as [nf^i,H] = taking values ei- 
ther 1 or 0, according to the whether site i is occupied or 
unoccupied by an f-electron. Due to the local conservation 



of /-electron number the Hamiltonian (1) may be written 
as, 

H = Y,h.,{uj)d+d,+EfJ2co, (2) 

where hij{uj) = —tij + (UuJi ~ ^)Sij. 

We set the hopping integral tij ~ 1 , for nearest neigh- 
boring sites i and / and = 0, otherwise. The eigenvalue 
spectrum of the Hamiltonian H, for a configuration uj of 
/-electrons is calculated by numerical diagonalization on 
a triangular lattice of finite size with periodic boundary 
conditions (PBC). The average value of physical quanti- 
ties is obtained by the classical Monte Carlo method using 
Metropolis algorithm. The details of the method can be 
found in our earlier papers [271130] . 

We work at half-filhng, i.e., Nj+Nd = N where Nt, Nd 
are the total number of /- and d-electrons and N = L {L = 
6) is the total number of sites. At a finite temperature T, 
the thermodynamic quantities are determined as averages 
over various configurations uj{Nf) with statistical weight 
P{u!{Nf)) given by 

-f!Fiuj{Nf)) 

Pi^m) = (3) 

where the corresponding free energy is given as. 
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Fig. 2. (color online) The d-electron energy spectrum for U=0. Fig. 3. Specific heat (C„) as a function of temperature calcu- 
lated at f/ = for different values of Ef. 
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(4) 



with ^ being the chemical potential. The /.t has been fixed 
to satisfy the condition n/+nrf=l. 
The partition function is 



N 



-tSEf. 



11(1 



-|3[\,{u^{Nf))-^,]^ 



(5) 



Let a thermodynamic quantity 'A' have value A{uj{Nf)) 
corresponding to the configuration u!{Nf) of A^/-electrons, 
then the ensemble average of 'A' at temperature T is ob- 
tained as 



(A) 



-PF{Lj{Nf)) 



(6) 



For example the ensemble average of number of f-electrons 
for given values of U, E/, T is obtained as 



{Nf) 



(7) 



and n/ = (uf) = 

The total internal energy E{uj) of the system correspond- 
ing to the configuration LolNf) is given as, 

EiuiNj)) = Y: + 1 + EfN.i.) (8) 



The specific heat (Ci,) is calculated from the fiuctuation- 
dissipation theorem (FDT). The specific heat is related 
through the FDT to the internal energy and defined as, 



(9) 



3 Results and discussion 

3.1 ISIoninteracting case (U=0) 

We firstly study the noninteracting case {U = 0). Using 
the method described above, we have found average value 

of f-electron occupation number per site n/ = (n/) = ^^^jp- 
as a function of Ef at different temperatures. Also we have 
calculated specific heat Cy as a function of temperature 
for different values of Ef. In Fig.l we have shown as 
a function oi Ef for different temperatures and also as a 
function of temperature for different values oi Ef (shown 
in inset of Fig.l). In Fig. 2 we show the d-electron energy 
spectrum corresponding to the configuration uj{Nf) with 
minimum internal energy. 

When Ef is changed upwards (starting from the bot- 
tom of conduction band), Uf is found to decrease suddenly 
at some particular positions oi Ef. For example at Ef=- 
4.0, Uf decreases abruptly from Uf ^ 0.98 to n/= 0.8, 
at £'/=-1.0, from 0.8 to 0.6 at £'/=0.0, from 0.6 to - 
0.48 and at £'/=2.0, from 0.48 to 0.03. If we look at 
the d-electron energy spectrum, we find d-levels at energy 
E=-4.0, -1.0, 0.0, 2.0; i.e. there are energy gaps [31] in d- 
elcctron spectrum in energy range [-4, -1], [-1, 0], [0, 2]. As 
nf+nd=l^ the chemical potential /i is pinned at Ef. So, as 
Ef is shifted upward, ^ shifts upward. When /i shifts from 
-4 to -1, there being no d-levels in this range, Ud remains 
constant and so should Uf. This is exactly what we find 
as shown in Fig.l. As ^ reaches the energy -1, where there 
are many d-levels available, these get occupied and n f de- 
creases abruptly. Now with the increase of temperature, 
d-states above chemical potential are also occupied frac- 
tionally. That is why Uf decreases a bit as temperature 
increases . However, in contrast to the case of bipartite 
lattice [32], temperature is unable to smear out the dis- 
continuous transition on the triangular lattice. Variation 
of with temperature at different Ef is shown in the 
inset of Fig.l. When Ef is near the bottom of the d-band, 
TLf is near unity at T — >■ and decreases slowly as T is 
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increased. At other positions oi Ef, nf is lesser than unity 
(and the system is in classical mixed- valent regime) and 
decreases slowly as T is increased. 

Shown in Fig. 3 is the specific heat Cy as a function of 
temperature for different positions oi Ef. A broad single 
peak is observed in specific heat curve at all values oi Ef 
except at Ef = —2.0. At Ef = —2.0 there are two peaks 
(one sharp peak followed by a broad peak). At low tem- 
perature we observe that variation in Cy deviates from the 
linear behavior of free electrons. The nature of the peaks 
in the specific heat curves could be understood if we com- 
pare these with the d-electron energy spectrum (shown in 
Fig. 2) for different temperatures. The d-electron spectrum 
does not depend on Nf, Ef and temperatures. The spec- 
trum is discrete and multi-gaped. Wide gap around Fermi 
energy in d-clcctron spectrum corresponds to broad peak 
and small gap corresponds to sharp peak in Cy. 



3.2 Finite interaction case {U ^ 0) 

Let us consider the case where Coulomb correlation U 
between d-and f- electrons is finite. Fig. 4(a), (b) and (c) 



show {nf — Ef) phase diagrams for U =1, 5 and 10 and for 
temperatures 0.01, 1, 2, 3 and 4 respectively. We observe 
that for 11=1, as temperature increases from 0.01 to 4, the 
valence transition (i.e. nj-transition) (i) occurs at higher 
values of Ef and (ii) becomes smoother. We observe that 
the transition width (the range oi Ef over which Uf goes 
from 1 to 0), increases as temperature increases from 0.01 
to 4. 

At U=5 and 10 the [rif — £'/) -transition at different 
temperatures shows behavior similar to U=l for Ef well 
below or well above the middle of the conduction band. 
When Ef \s near the middle of the conduction band, fur- 
ther increment in the temperature does not affect the 
n/-transition significantly. This phenomena could be ex- 
plained by the large gap in the many body states around 
the Fermi energy (Ep). Due to large gap at the Fermi en- 
ergy (see Fig. 7 (b) and (c) for U=5 and 10 respectively) 
even higher temperature is not sufficient to move the lo- 
calized f-clectrons to the higher d-statcs. 

Fig. 5(a), (b) and (c) show variation of nf with tem- 
perature for U =1, 5 and 10 and for different Ef positions. 
We observe that for all values oiU,nf approaches towards 
the value around 0.5 at high temperatures. At low tem- 
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Fig. 7. (color online) The d-electron energy spectrum for different values of temperature and Ef calcaluted at (a) U=l, (b) 
U=5 and (c) U=10. 



peratures Uf takes values m the range [0.5, 1.0] depending 
upon the values of U and Ef. Our results are in contrast 
to the results of P. Farkasovsky [32] obtained on bipartite 
lattice. We observed discontinuous {nf — Ef) transitions 
with temperature. We report that even higher tempera- 
ture is not sufficient to smear the discontinuous {uf ~ Ef) 
transition observed in ground state spinless FKM on tri- 
angular lattice PTIISO] . They observed smooth variation 
in as a function of temperature. The contrast between 
these two results are due to the choice of different lat- 
tices. In square (bipartite) lattice, energy spectrum of the 
d- electrons are symmetric about the middle of the con- 
duction band, while in triangular (non-bipartite) lattice 



the d-electron energy spectrum is asymmetric about the 
middle of the conduction band and large number of d- 
electron states lie in the small regime above the middle of 
the conduction band. 

In Fig. 6(a), (b) and (c) we have shown variation of Cy 
with temperature for U=l, 5 and 10 at different Ef po- 
sitions. In the small Coulomb correlation limit (say ^7=1) 
and Ef at the bottom of the conduction band [Ef = 
—3.0), a single broad peak in C„ is observed. Moving the 
f-electron energy level towards the middle of conduction 
band two peak pattern (one sharp followed by a broad 
peak) is observed. For Ef at center of conduction band 
a sharp peak followed by a low broad peak is observed. 
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These different peaks can be explained by the energy spec- 
trum of the d-electrons shown in Fig. 7(a). In this figure we 
have shown the d-electron energy spectrum corresponding 
to the configuration Lu{Nf) with minimum internal energy. 
The d-electron spectrum is multi-gaped and depends upon 
the temperature and Ef. For U=l, Ef=0 and for low tem- 
perature, a gapless spectrum is observed. At finite tem- 
perature, we observe a gap around chemical potential in 
d-electron spectrum. Moving Ef from middle to bottom 
of conduction band say Ef = —2, a large gap is observed 
at low temperature. The gap disappears at higher tem- 
peratures. The small gap in d-electron spectrum leads to 
a sharp peak while large gap corresponds to a broad peak 
in Cy. 

For U=5 and Ef at the bottom of conduction band a 
broad peak in Cy is observed. Width of peak is reduced 
by moving f-electron level towards middle of conduction 
band. Two broad peaks in C„ at Ef = —1.0 are observed. 
A single peak is observed at Ef=0. At U=10, a single peak 
is observed in C„ for all observed Ef values. The width 
of the peak in Cy decreases with moving Ef from bottom 
to middle of conduction band. The critical temperature 
(Tc) at which peak occurs in Cy, reduces as Ef moves 
from bottom to middle of conduction band. If we compare 
Tc for Ef=Q and at different U- values we observe that 
Tc goes on decreasing with increasing U, similar to the 
result observed by M. Maska [29]. These features are well 
explained on the basis of variation of d-electron spectrum 
with Ef and temperature shown in Fig. 7(b) and (c) for 
U=5 and 10 respectively. 

In conclusion, we have studied the spinless Falicov- 
Kimball model on a triangular lattice at finite tempera- 
ture. The nature of metal-insulator transition with tem- 
perature has been looked into in various parameter regime. 
The temperature dependence of specific heat (Cy) is also 
studied for different ranges of parameters like U, Ef. Dis- 
continuous metal-insulator transition is observed even at 
finite temperature. This is unlike the case of bipartite lat- 
tices where it has been reported that with temperature 
metal-insulator transition becomes second order in na- 
ture. But here we see that the temperature can not smear 
the discontinuous metal-insulator transition observed in 
ground state on triangular lattice. Single and double peak 
patterns are observed in Cy depending on the parameters 
U, Ef and temperature. It is proposed that the various 
features could be explained by many body spectrum of 
d— electrons. 

Acknowledgments. UKY acknowledges CSIR, India 
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